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ABSTRACT 



We present a numerical model describing a circularly symmetric gaseous disk 
r-| \ around the Be star x Ophiuchi. The model is constrained by long-baseline inter- 

♦^1 ferometric observations that are sensitive to the Ha Balmer line emission from 

O I the disk. For the first time our interferometric observations spatially resolve the 

'o^ \ inner region of the circumstellar disk around x Ophiuchi and we use these re- 

' suits to place a constraint on the physical extent of the Ha-emitting region. We 

demonstrate how this in turn results in very specific constraints on the parame- 
ters that describe the variation of the gas density as a function of radial distance 
from the central star. 



Subject headings: stars: emission-line, Be — stars: individual (x Oph) — tech- 



O ' niques: interferometric 

OO ■ 



1. Introduction 

Models of the circumstellar gas around some B-type stars have been constructed for 
a number of years. Because the recombination in the circumstellar gas produces a series 
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of emission lines, most commonly in the hydrogen Balmer series, such systems are usually 
denoted B-emission or Be stars for short. Model s describing the g eometry of this circum- 
stellar ma terial have varied from sp herical shells (JGehrz et al.lll974j ). to slabs with constan t 
thickness (JKastner fc Mazzalilll989l ). to disks with constant opening angles (IWaterslll986l ). 
to disks with an exponential density fall-off in the direction perpendicular to the plane of 
the disk (JMarlboroughl Il969l ) . Although the first models assumed a constant temperature 
throughout t he disk, more recent rnodels, such as the ones deve l oped by iMillar fc Marlborough 
( 1l998l . ll999l ). ICarciofi fc BjorkmanI (120061 ) . and lSigut fc Joned (120071 ) . obtain a self-consistent 
temperature throughout the disk by enforcing radiative equilibrium (namely, by balancing 
the microscopic heating and cooling rates operating within the disk). 

While the temperature structure of Be star disks is now the subject of sophisticated 
modeling, the density structure largely remains ad-hoc with simple, parameterized models. 
Perhaps the most popular is a power-law density decrease with radius in the equatorial 
plane coupled with the aforementioned exponential drop perpendicular to the disk. This 
vertical density structure is assumed to result from the gravitational equilibrium set by 
the balance between the gradient in the gas pressure and the vertical component of the 
star's gravitational acceleration. Despite the simplicity of such models, they have been quite 
successful in matching a wide range of observational diagnostics that are sensitive to the 
geometric distribution of the circumstellar gas. These inclu de the infrared ex cess exhibited 
by Be stars over normal B stars of the sa me spectral class (JDoughertyi Il994l ) . the observed 
continuum linear polarization of Be stars (IWood et al.lll997l). a nd the shapes of the emission 
lines observed in Be star spectra (see for example lHummelll2000l ) . All of these diagnostics rely 
only on spatially unresolved spectra or spectropolarimetry. However, even more powerful 
constraints on the geometry of the emitting regions in Be stars can be obtained by combining 
these spatially unresolved diagnostics with new optical interfero metric observations that 
directly resolve the structure of the emitting region on the sky (JQuirrenbach et al.l 119971 : 
Tvcner et allboOGbl : Icies et al.l[2007[ ). 



To demonstrate how specific constraints can be placed on a disk model using obser- 
vations that spatially resolve the source, we chose to concentrate on the Be star x Oph 
(B2Ve; HR 6118; HP 148184). T his relatively nearby star with a Hipparcos distance of 
150 ± 17 pc (jPerryman et al.l 119971 ) has a strong Ha emission and therefore is an ideal can- 
didate for probing the disk structure at these wavelengths. 
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Observations 



2.1. Interferometry 

We acquired interferometric observations of x Oph on four nights between 2006 June 11 
and 2006 June 18. All observations were made with the same instrumental configuration. 
We used 4 different telescopes, resulting in five unique baselines, with one baseline si gnal 
measured twice on two independent output beams (see Fig. 2 in iTycner et al.l l2006bl . for 
a schematic of the beam combiner). The baselines used and their corresponding lengths 
are: AC-AE (18.9 m), AC-AW (22.2 m), AW-W7 (29.5 m), AC-W7 (51.6 m), AE-W7 
(64.4 m), with the AC-W7 baseline measured on two output beams. The observations were 
obtained using two spectrographs, each recording simultaneously signal from 3 baselines at 16 
spectral channels in the wavelength range 560-870 nm. The resulting (m, f )-plane coverage 
of the interferometric observations obtained from the spectral channel containing the Ha 
emission line is shown in Figure [H 

The observations of x Oph were interleaved with observations of a nearby calibrator 
star C Oph (09V; HR 6175; HD 149757). The choice of the calibrator star was based on the 
relative proximity on the sky to the target and the spectral type, which was similar to that 
of the target star. Furthe r more, the angular diameter of C, Oph has been directly measured 
by iHanbury Brown et al.l (119741 ) using an intensity interferometer with sufficient accuracy 
to be competitive with other methods that estimate angular diameters based on theoretical 
or empirical models. Therefore, for the ca librator star, we adopted a u niform disk angular 
diameter of 0.50±0.05 mas as obtained by iHanbury Brown et al.l (119741 ). 



The data reductions followed standard NPOI procedures and will only briefly be de- 
scribed here. The raw data are processed to produce squared visibilities (V^^) averaged into 
1 s intervals. These data are then flagged to eliminate points with fringe tracking and point- 
ing problems (JHummel et al.lll998l . l2003l ). The flagged data set is then used to obtain 30 s 
averages (we refer to these averages as scans), which are bias corrected using data obtained 
off the fringe (these are also known as incoherent scans). Table [1] lists the individual nights 
and the number of scans acquired on each night. Lastly, the squared visibilities from the 
spectral channel c ontaining the Ha emission line are calibrated with respect to the contin- 



uum channels (see iTycner et al.l l2003l . l2006al . for a detailed description of the procedures) 



Our flnal data set (shown in Fig. [2]) consists of the ensemble of all the V"^ data from the Ha 
channel obtained on all four nights listed in Table [H The individual calibrated Ha V"^ values 
are also listed in Table [2J Because squared visibilities represent the normalized Fourier power 
of the source structure on the sky, the interferometric observations can be compared directly 
to the Fourier transform of the synthetic image produced by the model. 
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The observations of x Oph used in this study were not optimized for the purpose of 
cahbrating with respect to an external cahbrator. However, such a cahbration allows one 
to inspect the continuum channels for the signature of a resolved stellar disk or a binary 
companion. Using the scans obtained on C, Oph to obtain the instrumental system response, 
we found that large systematic residuals from scan to scan remained after the calibration, 
which we attribute to atmospheric variations on a timescale of minutes that are poorly 
sampled by the scans on the calibrator. For this reason, we can only conclude that the 
central star at the continuum channels is marginally resolved with a uniform disk diameter 
of < 1 mas. An interferometric signature of a resolved binary is less susceptible to systematic 
variations from scan to scan because of its distinct sinusoidal functional form that can be 
seen across different spectral channels. In fact, the NPOI has been used succe ssfully to detect 
binaries over a wide range of separations (see for example iHutter et al.ll2004l ). The search for 
a binary signature in observations of x Oph, however, did not yield any convincing evidence 
of such a signature in our interferometric observations that would correspond to a magnitude 
difference {Am) of ~ 3 or less at separations of few hundred milli-arcseconds (mas) or less. 
Although we cannot rule out a fainter companion, such a companion would not contribute 
significantly to the total emission from the system and therefore for the purpose of the 
analysis presented in this study, x Oph can be treated as a single star. It should be noted, 
however, that y Oph has been classified as a single-lined spectrosc opic binary by more 
than one study fJAbt fc Levylll978l : lHarmanedll987l : iLevato et al.lll987l ). althoug h there is no 
agreem ent on the pe riod of the binary , which was estimated to be 138.8 d by lAbt fc Levy 
fll978h and 34.1 d bv lHarmane3 Jl987f ). 



2.2. Spectroscopy 



To estimate the strength of the Ha-emission in x Oph at the time of interferometric 
observations we obtained a high resolution spectrum in the Ha region using a fiber-fed 
Echelle spectrograph at the Lowell Observatory's John S. Hall telescope. The spectrum 
was acquired on 2006 June 10, just one day before our one week long interferometric run. 
The spectroscopic data were processed using standard routin es developed specifically for the 
instrument used to acquire the observations (JHall et al.lll994J ). The final reduced spectrum in 
the Ha region reaches a resolving power of 10,000 and a signal-to-noise ratio of few hundred 
(see Fig. ED. 

In addition to the spectrum shown in Figure [3] we have also acquired spectroscopic 
observations of x Oph in 2006 May and 2006 September. The peak intensities and the 
equivalent widths (EWs) of the Ha emission line all agree to within 1% with the line profile 
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shown in Figure [31 Therefore, we assume that during our interferometric run the overall 
emission in the Ha line was stable and combining all of the interferometric observations into 
one data set is justified. 



3. Determining the Size of the Ha-emitting Region 



We followed the procedure described by iTycner et al.l (l2006bl ) to determine the angular 
extent of the Ha disk. We fitted a circularly symmetric Gaussian model to the interferomet- 
ric data and obtained a best-fit diameter (defined as full-width at half-maximum; FWHM) 
for the disk of 3.46 ± 0.07 mas, which at a distance of 150 pc corresponds to a diameter 
of 0.52 AU (112 Rq). The fit results in a reduced x^ value of 0.95, and therefore, there is 
no indication that the observational signature deviates from circular symmetry or a Gaus- 
sian shape for that matter (shown as dash-dotted line in Fig. [2]). In fact, this is the best 
observational evidence to date that a Gaussian model is not only the simplest mathematical 
description of the radial distribution of the Ha emission in a circumstellar disk, but it is also 
fully consistent with t he observed d i stribut ion. This is in addition to similar conclusions 
obtained previously by iTycner et al.l (l2006bl ) for two other Be stars, 7 Gas and (p Per. 



It is interesting to compare our resul ts based on the Gauss i an na odel fit with the simple 
estimator based on Ha EW obtained by iGrundstrom &: Giej (120061 ). With the stellar pa- 
rameters for X Oph listed in Tableland assuming that the disk is viewed at 20° (see ^15.21) . 
for an EW of —7.1 nm based on our spectroscopic observations IGrundstrom fc GiesI (120061 ) 
prediclo a disk radius of 9.2 R^,. This agrees well with the radius of ~ 9.8i?* we obtain based 
on the Gaussian fit (using half of our FWHM measure). 

Although characterizing the Ha emitting region with a Gaussian model can be a useful 
tool, especially when only a limited interferometric data set is available to constrain the 
disk characteristics, it is only a simple observational parameterization. However, in cases 
where the interferometric data set covers wide range of spatial frequencies it is desirable to 
extract more information about the physical characteristics of the disk, such as density and 
temperature distribution. For this purpose a much more sophisticated model is needed. 



^We have used the programs provided by the authors at'http://www.chara.gsu.edu/~gies/Idlpro/BeDisk.tar| 
to extrapolate for equivalent widths beyond the values considered in Fig. 1 of iGrundstrom fc Giea (2006;). 
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Disk Model 



The radiative equilibrium model used to represent the thermal structure of the gaseous 
circumstellar disk surr ounding y Oph has b een computed using the bedisk code, which is 
described in detail by ISigut &: Jonej (120071 ). This code incorporates many improvements 
over previous treatments, most notably the use of a solar chemical composition for the cir- 
cumstellar gas, which is an important ingredient for calculating heating and cooling rates. 
These models are ideally suited to compare with interferometric observations because model 
monochromatic, 2-dimensional images of the circumstellar disk at specific wavelengths can 
be computed. The bedisk code can also be used to compute the hydrogen line spectra and 
the overall spectral energy distribution. For the calculation of the hydrogen line profiles, the 
disk is assumed to be in pure Keplerian rotation. There is considerable evidence that Be 
star disks are indeed rotati onally supported based on the detailed analysis of line profiles 
( iHummel fc Vranckerul2000l ) and on the interpretation of V/R variations as o ne-armed den- 
sity waves in the disk, which are predicted for Keplerian disks (jOkazakil 120071 ). Recent work 
based on long-basel ine interferometry ha s also provided observational evidence in support of 
Keplerian rotation (JMeilland et al.ll2007l ). 



The disk density model of ISigut fc Joned (120071 ) requires only a few input stellar and 

disk parameters. It is assumed that the density in the equatorial plane of the disk is given 

by an R^"^ power-law in radial distance (where R is the distance from the star's rotation 

axis), and that the density distribution vertical to the equatorial plane (in the Z direction) is 

set by the hydrostatic equilibrium established by the gas pressure gradient and the vertical 

component of the star's gravitational acceleration. Such a model generally produces a thin 

disk in which Z/R <C 1, and in this case the form of the density distribution is particularly 

simple: 

/ p 

p{R, Z) = po 



R 



1 



where po is the density at the inner edge of the disk in the equatorial plane, n is the index 
in the radial power-law, and H is the scale height in the Z direction and is given by 



H 



'2i?3 



with the parameter ag of the form: 



"0 



GM, 



kTo ■ 



(2) 



(3) 



In these expressions, M^ and i?^ are the stellar mass and radius, respectively; Tq and Pq are 
the assumed (vertically isothermal) temperature and mean-molecular weight at the radial 
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distance R. This simple, analytical form for the density is possible because it is assumed 
that the vertical pressure scale height at each R can be represented by a single temperature 
Tq. Although this assumption might not be directly applic able to dense disk s where strong 
vertical temperature gradients in the disk might be present, ISigut et al.l (120071 ) examined the 
accuracy of this assumption in detail (in pure hydrogen disks) and found that in most cases, 
the differences in typical predicted diagnostics, such as the Ha profiles and IR excesses, 
are not large for models in consistent radiative and (vertical) hydrostatic equilibrium if the 
density-averaged disk temperature is used for the parameter Tq. The largest variations are 
found for the densest models (with largest po) as these models develop a cool, equatorial 
zone close to the star. Hence for this current work, we have adopted the simpler density 
structure given by equation [1] with the parameter Tq being constant for all values of R and 
Z. In the case of models constructed for x Oph we set Tq at 10,500 K based on typical 
density-weighted disk temperatures obtained for models in reasonable agreement with the 
observations of x Oph. 

With this model for the disk gas density the parameters po and n are free to be varied 
to match the interferometric visibilities and Ha line profile of x Oph. Before we can obtain 
a self-consistent solution for a disk model with specific po and n values, we need to adopt the 
parameters for the central star. The reported Tp ff values for y Oph in the li t erature range 



from the low of 18,000 K to a high of 29.600 K flGoraval ll984J: IWaterslll986l : iFremat et al. 



20051 : IZorec et al.ll2005l ). Furthermore, iGorayal (Il984j ) not only detected variations in visual 
magnitude of x Oph but also reported that such variations were accompanied by changes 
in the slope of the continuum, resulting in variations of derived T^^ values by more than 
4,000 K. Because we do not attempt to model the intrinsic variability of the source in this 
study, we adopt Tcs bas ed on the spectral type of y Oph, a main sequence B2 star, which 
based on tabulations of Ide Jager fc Nieuwenhuijzenl (119871 ) gives a T^q of 20,900 K. The 
stellar mass and radius are less critical parameters for t he di s k mod els and we adopted the 
average values based on the spectral type tabulated by ICoxl (120001 ) of 1O.9M0 and 5.7Rq. 
The model input parameters for the stellar component are listed in Table [31 

Figure m shows the thermal structure of a disk model with n = 2.5 and po = 2 x 10^^^ g 
cm~^ as a function of radial distance from the central star and distance from the equatorial 
plane. Although the computational grid extends out to 285 Rq (50 i?*) from the central 
star, the figure shows only the region out to 200 Rq. There is a conspicuous cooler region in 
the equatorial plane of the disk within the first ~ 50Rq from the central star. This cooler 
volume represents the region of the disk where the density and the corresponding optical 
depths are larger. The disk temperature reaches a minimum of ~ 6, 000 K in this volume. 
The density perpendicular to the equatorial plane falls off approximately exponentially, so 
regions surrounding this cool volume have significantly lower densities and are optically thin 



resulting in greater temperatures. These features in the disk thermal structure are typical of 
other models presented in the literature for Be stars with moderate to high di s k densities (for 
exara ple see IJones. Sigut &: Marlboroughl 12004 ICarciofi &: BjorkmanI l2006l : ISigut &: Jones 
20071 ). The overall average temperature of the disk over the computation grid is ~ 10, 400 K, 
which is consistent with our adopted temperature used in the vertical hydrostatic equihbrium 
calculations (recall Tq). It should also be noted that the temperature structure shown in 
Figure Hlis influenced by the incoming radiation from the central star and therefore ultimately 
the thermal structure will be affected by the choice of stellar parameters, such as T^g and 



5. Constraining the Model 
5.1. Synthetic H« Image 



A self-consistent temperature solution for the disk density model (like that shown in 
Fig. Hj) can be used to directly compute the specific intensity of radiation emitted perpen- 
dicular to the plane of the disk. As the hydrogen level populations are found naturally as 
part of the thermal solution, the Ha line emissivity and opacity, 

hu23 



vf 



and 



Xv 



hU23 

An 



Att 



{N2B; 



-iV^A 



3 ^32 (Pu , 



23 



N^B' 



32 J 



(4) 



(5) 



are known at each {R, Z) location in the disk. Here N2 and A^3 are the hydrogen number 
densities in levels 2 and 3 respectively, B23, -B32 and A32 are the Einstein probability coeffi- 
cients for Ha, and (j),^ is the Ha line profile. For the latter, we have adopted the routines of 
Barklem &: Piskunovl (120031 ) which include the contributions of (thermal) Doppler broadening 
and coUisional broadening (including the linear Stark effect). Turbulent (or microturbulent) 
broadening was not considered. 

To form the total opacity and emissivity at the frequency of Ha, we have added the 
continuous free-free and bound-free emissivity and opacity of hydrogen to equations H] and [51 
We have also included electron scattering as a coherent source of opacity but not as a 
source of emissivity. In a radiative equilibrium (or thermal-balance) calculation it is not 
usual to include continuum scattering as it leads to no net energy gain or loss by the gas. 
However, non-coh erent electron scattering is capable of broadening spectral lines and blurring 
continuum edges (JRybicki &: Hummerlll994l ). and this process could potentially effect the Ha 
EW measure. Although these effects are ignored in the current study, non-coherent electron 
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scattering will be included in future work. Given the opacity and emissivity, we solve the 
equation of radiative transfer at each radial distance R in the disk along a ray in the Z- 
direction, i.e., 



Ctl^yK, Zy 32 23 



dZ 



r^i'-xfh{R.Z). (6) 



The fo rmal solution was performed with the short-characteristics method of lOlson fc Kunasz 



(119871 ). This gives the specific intensity emerging perpendicular to the disk and is what an 



external observer would see if the disk was viewed at i = 0°. 

To obtain a profile that can be compared with observations, the kinematic broadening of 
the disk's (assumed) Keplerian rotation must be included. As we estimate that the inclination 
of the rotation axis of xOph is quite small {i ~ 20°, see below), we have assumed that the 
specific intensity at this small inclination is well approximated by the i = 0° intensities 
computed above. Given this, we simply shift the i = 0° spectrum of each projected area of 
the disk and star on the sky by its radial velocity. For areas on the st ellar surface, we have 



adopt ed the photospheric Ha profiles computed in LTE by the code of iBarklem fc Piskunov 



(120031 ). The radial velocity of each patch on the stellar surface followed from the assumption 
of a spherical star in solid body rotation with an equatorial velocity of 375 km s~^ (giving a 
V sin i of 144 km s~^; see below). Clearly this procedure is not appropriate for large inclination 
angles as the perpendicular rays through the disk would no longer be representative of the 
rays in the direction towards the external observer. However, direct comparison of the 
above procedure with a complete routine that solves the transfer equation along the inclined 
rays pointing at the observer shows that the simpler procedure gives good results for small 
inclinations and the approximate procedure is computationally much faster. 

We specify more than a hundred wavelength steps to cover the 15 nm wide spectral 
region centered at the Ha line, which corresponds to the spectral channel used to obtain the 
interferometric observations. The integrated net intensity over that 15 nm spectral window 
is shown in Figure [5] where clearly the total intensity is dominated by the bright central 
star. This is expected because the 15 nm wide spectral region is much wider than the Ha 
emission line (recall Fig. [3]) and therefore the Ha channel receives significant fraction of the 
light from the central star. In turn, this radial intensity distribution can be used to construct 
a circularly symmetric image that corresponds to the total intensity in the Ha channel. 



The relatively small v sin i value reported for x Oph of 144 km s by IZorec et al. 



(120051 ) and the lack of interferometric signature that would suggest large deviations from 
circular symmetry (recall §3]), both suggest that the inclination angle of the disk to the 
plane of the sky is small. Assuming that the star is rotating near its critical velocity of 



477 km s ^ (estimated based on its spectral type; see lPorterlll996l . and references therein), the 



inclination angle i must be ~ 18° or more for the star not to rotate above its critical velocity. 
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As rapid rotation is a known characteristic of Be stars, we adopt 20° for the inchnation 
angle of x Oph. Furthermore, because the axial ratio produced by the projection effect on 
geometrically thin and circularly symmetric disks scales as a cosine of the inclination angle, 
an inclination of 20° produces only a ~6% departure from circular symmetry (i.e., an axial 
ratio of 0.94). For circumstellar disks with non-negligible opening angles the effect would 
be even smaller. However, an axial ratio of 0.94 would be undetectable in our data because 
it affects the V'^ values at a level much smaller than the precision of our interferometric 
observations (compare the dash-dotted and dotted curves in Fig. [2]). Therefore, for the 
purpose of this study we will treat the synthetic images as circularly symmetric and will 
ignore the small projection effect on the images. 

The squared visibilities from the Ha channel measure the normalized Fourier power of 
the source structure on the sky, and therefore the Fourier transform of the model image and 
the interferometric observations can be directly compared. To accomplish this we first con- 
struct a circularly symmetric image corresponding to the radial distribution of the integrated 
intensity over the 15 nm wide spectral bandpass covered by the NPOI spectral channel con- 
taining the Ha emission line. Because the model intensities are calculated for radial region 
extending up to 50 i?^ from the central star, we require minimum dimensions of the model 
image of 100 i?^ x 100 i?^ to fully describe the output. However, we extend our synthetic 
image to 1000 Ri, x 1000 i?* dimensions (the outer regions are filled with zeros) to increase 
the sampling frequency of the model in the Fourier space. The image is sampled every 0.2 i?^ 
and this is sufficient to avoid any aliasing problems at the high spatial frequencies covered 
by the observations. The region of interest of the synthetic image is shown in Figure Ei 
The 2-D Fourier transform of this circularly symmetric image results in a function that also 
has circular symmetry. The normalized Fourier power of this transform can then be plotted 
as a function of the radial spatial frequency (shown as a solid line in Fig. [2]). We use the 
standard x^ statistic to assess the goodness-of-fit of the model to the actual data obtained 
on the source. 



5.2. The Parameter Space 

We explore the parametric region of of the disk density model (eq. [1]) by computing over 
500 parameter pairs of n and po- We covered the range of 1.8 - 5.3 in n and adjusted po from 
a low of 1.0 X 10'^^ g cm~^ to a high of 8.0 x 10~^ g cm~^. The majority of the solutions 
were completely inconsistent with the interferometric observations, and those solutions that 
yielded reduced x^ value of 7 or less were concentrated in a very well defined 'valley' in the 
X^ space (see Fig. [7]). 
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The model that results in the lowest reduced x^ value of 1.17 corresponds to n of 2.5 and 
Po of 2.0 X 10"^^ g cm~^. The thermal structure (Fig. H]), Ha model image (Fig. [6]), and the 
model curve shown in Figure [2] all correspond to this best-fit model. Furthermore, assuming 
that the disk is viewed at an inclination angle of 20° (as discussed in § 15.11) . we obtain an 
excellent agreement with the observed Ha line profile obtained through spectroscopy (see 
Fig. [3]). In fact, for the disk model defined by our best-fit to interferometric data, the Ha 
profile can only be reproduced assuming i of 20°. For any smaller inclination the model 
profile is simply too sharp (due to the decreased rotational broadening) and for larger values 
of i the line becomes too broad and weaker than the actual observed profile. 

In addition to the best-fit solution found at n of 2.5 and po of 2.0 x 10^^^ g cm~^ it is 
clearly evident from Figure [7] that there exists a range of parameter values in the x^ space 
that produce disk models that fit the observational data acceptably well. However, these 
solutions appear to be confined to a very narrow range in the x^ space. We attribute this 
range of solutions to the changing total Ha flux generated by a model, where most of the 
solutions in the upper-left part of the figure produce too much Ha flux with respect to the 
central star and those in the lower-right part of the figure do not produce enough Ha flux. 
These effects can be related to different disk thermal structures and different densities, which 
in turn will affect the apparent sizes of the Ha-emitting regions. A linear least-squares fit to 
solutions with reduced x^ < 7 in the log po versus n plane yields a slope of 1.18±0.03 (and an 
intercept of — 13.6±0.1), where the uncertainty does not account for the effects of gridding or 
the choice of the reduced x^ cutoff value. Nevertheless, this indicates that at least for a range 
of n values between 2.5 and 4.0 the disk solutions that produce interferometric signatures 
closest to the observational data fall along a relation wh ere log pn is directly proportional to 



n. This is in agreement with similar conclusions made by lGies et al.l (120071 ) who approximate 
the size of the disk on the sky using a boundary between the optically thick and thin disk 
regions, which then leads to the expectation that good fits fall along a relation where log po 
is proportional to n. 



6. Discussion 



It is useful to compare our best- fit disk model for x Oph to previous results of [Waters 



( 1l986l ) who modeled the IRAS IR excesses of a large sample of Be stars with disk models 
with fixed opening angles of 15° and where the density varied only with the distance from 
the stellar center, r, in the form of 

p{r) = p'o{r/R.)-r (7) 
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Assuming an isothermal disk with a temperature of 18,000 K, a temperature that is con- 
siderab ly larger than the ~10,400 K density-weighted average temperature of our best-fit 
model, IWaterd (Il986l ) found a radial power-law index n of 2.4 and a density at the base 
of the disk of logpo = —11-4. These results compare well to our results oi n = 2.5 and 
logpo = —10.7, even though the disk-density models used in the two studies are not the 
same. The smaller base disk density found by IWaterd (119861 ) is likely the result of the fact 
that his density distribution is constant along radial arcs over the entire opening angle of 
the disk which places more material in the disk as compared to our models. 



Porterl (119991 ) modeled the circumstellar regions of Be stars as isothermal viscous disks 



and predicted n of 7/2. Viscous models that drop the assumption of an is othernial disk 
can find a range of n values as the disk fills or empties (I Jones et al.l 120081) . iPorterl (119991 ) 



also analyzed the IR excess of x Oph using the data of IWaterj (119861 ) . but with a disk 
model den sity ess e ntially equivalent to ours. In the case of an isothermal disk (also set to 
18,000 K), iPorterl (119991 ) found an index n of 2.2 with logpo = —11-2; in the case of a disk 
in which t he dens i ty-ave raged temperature was allowed to vary with radial distance as a 
power-law, IPorterl (1l999l ) obtained n = 1.9 and log(po) = —11-2. Porter indicates that the 



non-isothermal disk formally fits the IR data best, but that the improvement is not large. 

In addition to the comparison between model and observed Ha interferometric visibilities 
and the Ha line profile, it is also instructive to compare the predicted IR excess of our best-fit 
model to the IR observations reported in the literature. Figure [H] illustrates the com parison 
between the SED of our best-fit model and the visual and IR observations from Waters 



(119861 ). The overall agreement is quite reasonable, especially considering that the model was 
not fit to the SED data, although the observations in the 1 — 10 /xm region fall be low the 
model. We should also point out that variability between the observations used by [Waters 
(119861 ) and o ur 2006 Ha observation s (interferometric and spectroscopic) cannot be ruled out. 
For example, iBanerjee et al.l (120001 ) obtained the Ha profile for x Oph in 1998 and reported 
an equivalent width of approximately —3.6 nm, as compared to the —7.1 ± 0.2 nm found 
in the present work for the 2006 epoch (recall Fig. [3]). Hence, the Ha emission of xOph 
has approximately doubled in 8 year s, and assumi ng that the Ha emission was similarly 
weaker at the time of observations of IWaterd (119861 ) we can then expect our best-fit model 
to over-predict IR excess as compared to Waters' results. For this reason, we feel that our 
model cannot be improved using the IR observations that were acquired using the IRAS 
satellite in the early 1980's. 
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7. Summary and Future Work 



In this paper we presented the study of the circumstellar disk of the Be star x Oph, 
based on the approach of combining high-spatial resolution interferometric observations with 
numerical disk models requiring only a few input parameters. Using the interferometric data 
we were also able to determine that the Ha-emitting disk can be fit by a circ ular symmetric 



Gaus sian with FWHM diameter of 3.46±0.07 mas. Using the code bedisk (jSigut fc Jones 



20071 ) we created a grid of models with a range of inner edge densities (po) and equatorial 
plane density distribution fall-off rates (n), while keeping the mass, radius and effective 
temperature of the central star fixed. These models were used to create synthetic Ha images 
by adding the emission line and continuum flux, including the contribution from the central 
star, in a 15 nm spectral region around Ha, similar to that of the Ha interferometric channel. 
These images were then Fourier transformed and compared directly to the interferometric 
observations, resulting in the best-fit model with n = 2.5 and po = 2.0 x 10~^^ g cm~^. The 
best-fit disk model was also used to calculate a synthetic Ha line profile and spectral energy 
distribution of this system over 0.1-100 pm, which show good agreement with the observed 
line profile and the photometric observations published in the literature. 

The technique demonstrated here, of combining models with interferometric observa- 
tions, demonstrates a new and independent method for obtaining properties of circumstellar 
disks of Be stars. This is possible because the density distribution in the disk is strongly 
affected by the value of the disk density at the stellar surface, po, and to a lesser extent by n. 
Variations in po and n result in dramatic changes in the thermal structure of the disk. Ul- 
timately, the density distribution and the resulting thermal structure directly affect the Ha 
flux predicted by the theoretical models, which in turn can be constrained by observations 
that spatially resolve the region. In the future, we plan to investigate more sophisticated 
forms of the model density distri bution within Be star d isks, such as rotating density pertur- 



bation models (see for example IWisniewski et al.l 120071 ). and density distributions that are 



consistent with hydrodynamic simulations. 

The Navy Prototype Optical Interferometer is a joint project of the Naval Research 
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Table 1: NPOI Observations of x Qph 

UT Date 7^ of Scans 7^ of Baselines 

2006 June 11 5 5 

2006 June 13 4 5 

2006 June 17 6 5 

2006 June 18 7 5 



Table 2: Calibrated Ha Squared Visibilities of x Oph 



Julian Date 


Spatial Frequency u 


Spatial Frequency v 






(JD - 2,450,000) 


(10^ 


cycles/radian) 


(10^ 


cycles/radian) 


VL 


Baseline^ 


3897.785 






27.003 






-7.122 


0.749 ± 0.048 


AC-AE 


3897.785 






64.880 






23.478 


0.534 ± 0.032 


W7-AC 


3897.785 






91.883 






16.357 


0.502 ± 0.023 


W7-AE 


3897.785 






-29.723 






-8.345 


0.688 ± 0.061 


AC-AW 


3897.785 






64.989 






23.518 


0.519 ± 0.027 


W7-AC 


3897.785 






35.266 






15.173 


0.680 ± 0.030 


W7-AW 


3897.801 






27.324 






-7.970 


0.696 ± 0.042 


AC-AE 


3897.801 






61.798 






21.501 


0.548 ± 0.045 


W7-AC 


3897.801 






89.123 






13.531 


0.520 ± 0.026 


W7-AE 


3897.801 






-28.537 






-7.436 


0.724 ± 0.043 


AC-AW 


3897.801 






61.902 






21.537 


0.511 ± 0.027 


W7-AC 


3897.801 






33.365 






14.101 


0.681 ± 0.037 


W7-AW 



Note. — Table[2]is published in its entirety in the electronic edition of the Astrophysical Journal. 
■f There are 6 baselines per scan (same JD) with the W7-AC measured at two output beams with 
slightly different wavelength scales and thus resulting in slightly different u and v values for the same 
baseline length. 
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Table 3: Model Stellar Parameters 

Parameter Symbol Value Reference 

Mass (Mq) W^ HTg Cox ('2000) 

Radius (i?©) R^ 5.7 Cox QOOO) 

Effective temperature (K) .... T^s 20,900 de Jaeer &: Nieuwenhuiizen Q 

Luminosity (L©) L^ 5.6 x 10^ calculated''' 

Surface gravity ^ogg 4.0 calculated* 

Note. — ^ L^, = AnRlaT^g where a is the Stefan-Bohzmann constant. * \ogg = log(G'M^i?~^) 
where G is the gravitational constant. 
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Fig. 1. — The (M,f )-plane coverage of the interferometric observations of x Oph obtained 
in the Ha channel with five unique basehnes from four nights of observations (circles). The 
dotted-hnes indicate possible coverage from meridian to 6 hr east and the solid-lines indicate 
possible coverage from the meridian to 6 hr west. 
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Fig. 2. — 5?5terfefi5liMEfe7d^Mt<febtaiae^M(m)the spectral channel containing the Ha emission 
line from the circumstellar disk of x Oph. The 0.52 AU diameter (FWHM) circumstellar 
disk is fully resolved, whereas the central star is assumed to have an uniform disk diameter 
of 0.35 mas [dashed line). The best-fit model [solid line) obtained by taking the Fourier 
transform of the synthetic image from Figure [6] is shown along with a circularly symmetric 
Gaussian model with FWHM diameter of 3.46 mas [dash-dotted line). The effect of an axial 
ratio of 0.94 on the Gaussian model is also shown [dotted line). 
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Fig. 3. — The Ha emission line of x Oph observed on 2006 June 10 {solid-line) and the 
synthetic Ha profile calculated from the best-fit disk model (dashed-line) . The model profile 
has been obtained for an inclination angle i of 20° and has been broadened by a Gaussian 
kernel to match the resolving power {R = 10,000) of the spectroscopic observations. The 
observed Ha profile has an EW of —7.1 ± 0.2 nm and the theoretical profile has EW of 
—6.6 nm. 
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Fig. 4. — The thermal structure of the disk model of x Oph obtained using n = 2.5, 
Po = 2.0 X 10~^^ g cm~^ and the stellar parameters listed in Table [31 The thermal structure 
is shown as a function of radial distance {R) from the center of the star and the vertical 
distance {Z) above and below the equatorial plane. 
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Fig. 5. — Integrated model intensity for the direction normal to the plane of the disk in a 
15 nm wide spectral region centered at Ha [solid line) for the best-fit model of n = 2.5 and 
Po = 2.0 X 10^^^ g cm~^. The Gaussian radial intensity distribution that reproduces the 
interferometric signature in Fig. [2] is also shown [dashed line). 
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Fig. 6. — Synthetic image of the best-fit disk model viewed pole on. Only intensities in a 
15 nm spectral window centered at Ha emission line were used to construct this synthetic 
image. The surface brightness of the central star is derived using the stellar atmosphere 
model integrated over the same 15 nm spectral window as used to calculate the net intensity 
from the disk. 
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Fig. 7. — Contour plot showing the reduced x^ value as a function of po (g cm"^) and 
n. Only a region of reduced x^ < 7 is shown and the best-fit model with reduced x^ of 
1.17 is marked with a plus sign. Some of the gaps in the linear trend are caused by the grid 
spacing and the limitations of the contour plotting routine dealing with unequal grid spacing. 
The approximate boundary of the region of parameter space sampled by our models is also 
indicated [dotted lines), along with a linear fit {dashed line) to the models with reduced 
X^ < 7 in the log po and n plane (see § 15.21) . 
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Fig. 8. — The spectral energy distribution (SED) of our best-fit model corresponding to 
n = 2.5 and po = 2.0 x 10^^^ g cm^^ (solid-line), which consists of the sum of the stellar 
contribution (dotted-line) and IR excess produced by the circumstellar disk. The mo del SEP 
is co mpared to visual, near-IR and far-IR observations reported in literature (see IWaters 
19861 . and references therein). 



